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I ' Major advances have been made in the field of gamma-ray bursts in the last two years. 

The successful discovery of X-ray, optical and radio afterglows has made possible the iden- 
tification of host galaxies at cosmological distances. The energy release inferred in these 
outbursts place them among the most energetic and violent events in the Universe. They 
are thought to be the outcome of a cataclysmic stellar collapse or compact stellar merger, 
leading to a relativistically expanding fireball, in which particles are accelerated at shocks 
and produce nonthermal radiation. The substantial agreement between observations and the 
theoretical predictions of the standard fireball shock model provide confirmation of the basic 
I | aspects of this scenario. New issues being raised by the most recent observations concern the 

amount and the nature of the collimation of the outflow and its implications for the energet- 

■ ics, the production of prompt bright flashes at wavelenghts much longer than gamma-rays, 
^.j. ' the time structure of the afterglow, its dependence on the central engine or progenitor system 

, behavior, and the role of the environment on the evolution of the afterglow. 

(N ; 

' §1. Introduction 

ON '■ 

Gamma-ray bursts (GRB) are brief pulses of gamma-rays which pierce, for a 
Qh' brief period of tens of seconds, an otherwise pitch-black gamma-ray sky. They are 

detected about once a day, and while they are on, they outshine everything else in the 
gamma-ray sky, including the Sun. In fact, they are the brightest explosions in the 
Universe. Until a few years ago, they left no trace at any other wavelengths besides 
gamma-rays, until in early 1997 the Italian-Dutch satellite Beppo-SAX suceeded in 
• i-H ! providing accurate X-ray measurements that allowed, after a delay of some hours 

to process the position. This made possible their follow-up with large ground-based 
optical and radio telescopes, which paved the way for the measurement of redshift 
distances, the identification of candidate host galaxies, and the confirmation that 
they were at cosmological distances, comparable to that of the most distant galaxies 
and quasars ever measured. Since even at those tremendous distances they appear 
so bright, their energy output needs to be stupendous. It is comparable to burning 
up the entire mass-energy of the Sun in a few tens of seconds, or to emit over that 
same period of time as much energy as our Milky Way does in a hundred years. The 
current interpretation of this radiation is that this large amount of energy, released 
in a very short time in a very small region, expands in highly relativistic outflow, 
which undergoes both internal shocks producing gamma-rays, and later develops a 
blast wave and reverse shock, as it is decelerated by interaction with the external 
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medium. The remarkable thing about this model is that it successfully predicted the 
expected observational properties of the afterglows which such bursts leave at X-ray, 
optical and radio wavelengths, which can be observed over months. This fireball 
shock model and the blast wave model of the ensuing afterglow have become the 
leading paradigms for the current understanding of GRB. 

GRBs were first discovered in 1974 by the Vela military satellites, monitoring 
for nuclear explosions in violation of the Nuclear Test Ban Treaty. When these 
mysterious gamma-ray flashes were first detected, which did not come from the 
Earth's direction, the first suspicion (quickly abandoned) was that they might be the 
product of an advanced extraterrestrial civilization. Soon, however, it was realized 
that this was a new and extremely puzzling cosmic phenomenon. For the next 25 
years, only these brief gamma-ray flashes were observed, which vanished too soon to 
get an accurate angular position and hence left no trace, or so it seemed. Gamma- 
rays are notoriously hard to focus, so no gamma-ray "images" exist to this day: 
they are just pin-points of gamma-ray light. This led to a huge interest and to 
numerous conferences and publications on the subject, as well as to a proliferation 
of theories. In one famous review article at the 1975 Texas Symposium on Relativistic 
Astrophysics, no fewer than 100 different possible theoretical models of GRB were 
listed, most of which could not be ruled out by the observations then available! 

A major advance occurred in 1992 with the launch of the Compton Gamma- 
Ray Observatory, whose superb results are summarized in a review by Fishman & 
Meegan 29 ) . The all-sky survey from the BATSE instrument showed that bursts were 
isotropically distributed, strongly suggesting a cosmological, or possibly an extended 
galactic halo distribution, with essentially zero dipole and quadrupole components. 
The spectra are non-thermal, typically fitted in the MeV range by broken power- 
laws whose energy per decade vF v peak is in the range 50-500 KeV 3 \ the power 
law sometimes extending to GeV energies 45 ^ . GRB appeared to leave no detectable 
traces at other wavelengths, except in some cases briefly in X-rays. The gamma-ray 



Fig. 1. Time Profile of a Typical 
Gamma Ray Burst. The y-axis is 
the photon count rate in the 0.05-0.5 
MeV energy, the x-axis is the time in 
seconds after the trigger. Both be- 
fore and after the burst, no gamma- 
rays are detectable from the direc- 
tion of the burst (Fishman & Mee- 
gan, 1995). 
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durations range from 10 -3 s to about 10 3 s, with a roughly bimodal distribution of 
long bursts of i&>,2 s and short bursts of tb<2s 55 \ and substructure sometimes down 
to milliseconds. The gamma-ray light curves range from smooth, fast-rise and quasi- 
exponential decay (FREDs), through curves with several peaks, to highly variable 
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curves with many peaks 29 )> 56 ). The pulse distribution is complex 103 )' 82 ), and the 
time histories can provide clues for the geometry of the emitting regions 27 )> 28 ). 

At cosmological distances the observed GRB fluxes imply enormous energies, 
and from causality they must arise in a small volume in a very short time, so an 
— 7 fireball must form 88 )' 39 )' 131 ), which would expand relativistically. The main 
difficulty with this was that a smoothly expanding fireball would convert most of its 
energy into kinetic energy of accelerated baryons (rather than photons), and would 
produce a quasi-thermal spectrum, while the typical timescales would not explain 
events much longer than milliseconds. This problem was solved with the introduction 
of the "fireball shock model" m )> 66 ) ; based on the realization that shocks are likely 
to arise, at the latest when the fireball runs into an external medium, which would 
occur after the fireball is optically thin and would reconvert the kinetic energy into 
nonthermal radiation. The complicated light curves can be understood in terms of 
internal shocks 112 ) in the outflow itself, before it runs into the external medium, 
caused by velocity variations in the outflow from the source, 

The next major advance came in early 1997 when the Italian-Dutch satellite 
Beppo-SAX succeeded in providing accurate X-ray measurements which, after a 
delay of 4-6 hours for processing, led to positions 17 ) , allowing follow-ups at optical 
and other wavelengths, e.g. 141 ). This paved the way for the measurement of redshift 
distances, the identification of candidate host galaxies, and the confirmation that 
they were indeed at cosmological distances 78 )' 23 )' 57 ). The detection of other GRB 
afterglows followed in rapid succession, sometimes extending to radio 30 )' 31 ) and over 
timescales of many months 142 ) , and in a number of cases resulted in the identification 
of candidate host galaxies, e.g. 122 )' 11 )' 84 ), etc. The study of afterglows has provided 
strong confirmation for the generic fireball shock model of GRB. This model in fact 
led to a correct prediction 70 ), in advance of the observations, of the quantitative 
nature of afterglows at wavelengths longer than 7-rays, which were in substantial 
agreement with the data 143 )' 135 )' 146 )> 114 )' 152 ). 

A major issue raised by the large redshifts, e.g. 57 )' 58 ), is that the measured 7-ray 
fluences imply a total energy of order 10 54 (^ 7 /4vr) ergs, where Af2y is the solid angle 
into which the gamma-rays are beamed. A beamed jet would clearly alleviate the 
energy requirements, but it is only recently that tentative evidence has been reported 
for evidence of a jet 58 )' 33 )' 16 ). Whether a jet is present or not, such energies are 
possible 71 ) in the context of compact mergers involving neutron star-neutron star 
(NS-NS ) or black hole-neutron star (BH-NS ) binaries, or in hypernova/collapsar 
models involving a massive stellar progenitor 91 )' 108 ). In both cases, one is led to 
rely on MHD extraction of the spin energy of a disrupted torus and/or a central 
fast spinning BH, which can power a relativistic fireball resulting in the observed 
radiation. 

While it is at present unclear whether there is a single or several classes of GRB 
progenitors, there is a general consensus that they would all lead to the generic 
fireball shock scenario mentioned above. Two recent trends in the observations 
are the increasing hints pointing at a massive progenitor, at least for those bursts 
where an afterglow is detected 92 ), and the possibility that, at least in some bursts, 
a peculiar supernova may may be involved, e.g. 15 °) . Much of the current effort is 
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dedicated to understanding the different possible progenitors more specifically, and 
trying to determine how this affects the details of the observable burst and afterglow. 
Some of these progenitors, clearly, could be amenable also to observation with the 
current generation of gravitational wave detectors, and the numerical simulations of 
the collapse or merger share many interesting aspects with calculations motivated 
by gravitational wave searches. 

§2. Progenitors and Black Hole Systems 

In the last few years it has become apparent that most plausible GRB progenitors 
would be expected to lead to a central black hole (BH) and a temporary debris torus 
around it. This includes both the currently popular massive progenitor systems, 
such as hypernova or collapsars (including failed supernova lb [SNe lb], single or 
binary Wolf-Rayet [WR] collapse, etc. 153 )' 91 )) ")' 153 )' 71 )' 91 )' 108 ), NS-NS or NS-BH 
mergers 90 )' 71 ), Helium core - black hole [He/BH] or white dwarf - black hole [WD- 
BH] mergers 108 -*' 155 ), accretion- induced collapse 144 )' 109 ), etc.. An important point 
is that the overall energetics from these various progenitors do not differ by more 
than about one order of magnitude 75 ). Another possibility is massive black holes (~ 
10 3 — 1O 5 M0) in the halos of galaxies. Some related models involve a compact binary 
or a temporarily rotationally stabilized neutron star, perhaps with a superstrong 
field, e.g. 139 )' 137 )' 144 )' 133 ), which ultimately also should lead to a BH plus debris 
torus. 

The two large reservoirs of energy available in these systems are the binding 
energy of the orbiting debris, and the spin energy of the black hole 71 ). The first can 
provide up to 42% of the rest mass energy of the disk, for a maximally rotating black 
hole, while the second can provide up to 29% of the rest mass of the black hole itself. 
However, as in the related AGN case, the question is how to to extract this energy 
and how to deposit it at larger distances, where it produces optically thin radiation. 

A classical energy extraction mechanisms is the vv — > e + e~ process 25 ), which can 
tap the thermal energy of the torus produced by viscous dissipation. To be efficient, 
the neutrinos must escape before being advected into the hole; on the other hand, 
the efficiency of conversion into pairs (which scales with the square of the neutrino 
density) is low if the neutrino production is too gradual. Typical estimates suggest 
a fireball of <10 51 erg 119 )' 34 )' 62 ), except perhaps in the collapsar case where 108 ) 
estimate 10 52 ' 3 ergs for optimum parameters. If the fireball is collimated into a 
solid angle fij then of course the apparent "isotropized" energy would be larger by a 
factor {At:/ fij) , but unless fij is ^10~ 3 — 10 -4 this would fail to satisfy the apparent 
isotropized energy of 4 x 10 54 ergs deduced for GRB 990123 58 ) . 

A more efficient mechanism for extracting and transporting the energy from 
the vicinity of the BH may be through magnetic torques. Strong magnetic fields 
may be generated by the differential rotation in the torus 90 )' 81 )' 71 )' 50 ), in any of 
the above progenitors. In NS-NS mergers, even before the BH forms the shrinking 
binary may lead to winding up of the fields and dissipation in the last stages before 
the merger 64 )' 143 ). An even larger energy source is provided by the spin of the hole 
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itself, especially if formed from a coalescing compact binary, which is likelier to be 
rotating close to the maximal rate. Since the BH is more massive, it could contain 
more energy than the torus. The energy extractable in principle through MHD 
coupling to the rotation of the hole by the B-Z (Blandford & Znajek 7 )) effect could 
then be even larger than that contained in the orbiting debris 71 )' 91 ) . Collectively, 
any such MHD outflows have been referred to as Poynting jets. 

The various progenitors differ only slightly in the mass of the BH; they can 
differ more on the debris torus mass, can differ even more markedly in the amount 
of rotational energy contained in the BH. Strong magnetic fields, of order 10 15 G, 
are needed needed to carry away the rotational or gravitational energy in a time 
scale of tens of seconds 139 )' 137 \ which may be generated on such timescales by 
a convective dynamo mechanism, the conditions for which are satisfied in freshly 
collapsed neutron stars or neutron star tori 24 )' 53 -*. If the magnetic fields do not 
thread the BH, then a Poynting outflow can at most carry the gravitational binding 
energy of the torus. For a maximally rotating BH this is 42%, and for a slow-rotating 
BH this is 6% of the torus rest mass, multiplied by an additional efficiency factor 
e for converting MHD energy into radiation. The torus or disk mass in a NS-NS 
merger is 121) M d ~ 1CT 1 - 1O~ 2 M , and for a NS-BH, a He-BH, WD-BH merger or 
a binary WR collapse it may be estimated at 91) ' 34) M d ~ 1M . In the HeWD-BH 
merger and WR collapse the mass of the disk is uncertain due to lack of calculations 
on continued accretion from the envelope, so 1M Q is just a rough estimate. The 
maximum torus-based MHD energy extraction is then 

E max ,t ~ 1 - 10 x W 53 e(M d /M Q ) ergs . (2-1) 

If the magnetic fields in the torus thread the BH, the spin energy of the BH can 
in principle be extracted via the 7 ) (B-Z) mechanism ( 71 )). The extractable energy 
is 

E max ,bh ~ (f(a)/0.29)eM bh c 2 ~ 5 x 10 53 f(a)e(M bh /M Q ) ergs. (2-2) 

where f(a) = 1 - ([1 + VY^?]/2) 1 / 2 < 0.29 is the rotational efficiency factor, and 
a = Jc/GM 2 is the rotation parameter, which equals 1 for a maximally rotating 
black hole. The f(a) rotational factor is is small unless a is close to 1, where it rises 
sharply to its maximum value /(l) = 0.29, so the main requirement is a rapidly 
rotating black hole, a>,0.5. For a massive progenitor it is then imperative that the 
collapsed rotate fast, and hence a system with a compact binary companion seems 
called for. A binary fast-rotating WR scenario might do this, which does not differ 
much in its final details from the He-BH merger, depending on what fraction of the 
He core gets accreted along the rotation axis as opposed to along the equator 34 ) . For 
a fast rotating BH of 2.5 — 3Mq threaded by the magnetic field, the maximal energy 
carried out by the jet is then similar or somewhat larger than in the NS-NS case. 
Rapid rotation is also guaranteed in an NS-NS merger, since the radius is close to 
that of a black hole and the final orbital spin period is close to the required maximal 
spin rotation period. Some scenarios less likely to produce a fast rotating BH are 
the NS-BH merger (where the rotation parameter could be limited to a < M ns jM^ 
unless the BH is already fast-rotating), and the (single) failed SNe lb, where the 
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last material to fall in would have maximum angular momentum, but the material 
that was initially close to the hole has less angular momentum. Some examples are 
discussed in the accompanying article by Ruffert 12 °) . The main point that is worth 
stressing is that the total energetics differ between the various models at most by a 
factor 20 for a Poynting (MHD) jet powered by the torus binding energy, and by a 
factor of a few for Poynting jets powered by the BH spin energy, depending on the 
rotation parameter. 

The major difference between the various models is expected to be in the loca- 
tion where the burst occurs relative to the host galaxy (see §7). They are also likely 
to differ substantially in the efficiency of producing a directly observable relativistic 
outflow, as well as in the amount of collimation of the jet they produce. The con- 
ditions for the efficient escape of a high-/ 1 jet are less propitious if the "engine" is 
surrounded by an extensive envelope. In this case the jet has to "punch through" the 
envelope, and its ability to do so may be crucially dependent on the level of viscosity 
achieved in the debris torus (e.g. 62 )), higher viscosities leading to more powerful 
jets. The simulations, so far, are nonrelativistic and one can only infer that high 
enough viscosities can lead to jets capable of punching though a massive (several 
M Q ) envelope. This is facilitated, of course, if the envelope is fast-rotating, as in 
this case there is a centrifugally induced column density minimum along the spin 
axis, which might be small enough to allow punch-through to occur. If they do, a 
very tightly collimated beam may arise. "Cleaner" environments, such as NS-BH or 
NS-NS merger, or rotational support loss/accretion induced collapse to BH would 
have much less material to be pushed out of the way by a jet. In these, however, 
there is no natural choke to collimate a jet, which might therefore be somewhat wider 
than in a massive progenitor case. 

§3. The Generic Fireball Shock Scenario 

No matter what the progenitor is, the ultimate result must almost unavoidably 
be an e 1 * 1 , 7 fireball, which initially will be very optically thick. The initial dimensions 
must be of order ct var <10 7 cm, since variability timescales are t var >^l0~ s s. Most 
of the spectral energy is observed above 0.5 MeV, hence the 77 — > mean free 
path is very short. Due to the highly super-Eddington luminosity, this fireball must 
expand. Since many bursts show spectra extending above 1 GeV, the flow must be 
able to avoid degrading these via photon-photon interactions to energies below the 
threshold m e c 2 = 0.511 MeV 41 ). To avoid this, it seems inescapable that the flow 
must be expanding with a very high Lorentz factor F, since then the relative angle at 
which the photons collide is less than F _1 and the threshold for the pair production 
is diminished. This condition is 

^10 2 [(e 7 /10GeV)(e t /MeV)] 1 / 2 , (3-1) 

in order for photons e 7 to escape annihilation against target photons of energy et ~ 
1 Mev 63 )' 41 ). I.e., a relativistically expanding fireball is expected. From general 
considerations 67 ), an outflow arising from an initial energy E Q imparted to a mass 
M D << Eq/c 2 within a radius 77 will lead to an expansion where initially the bulk 
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Lorentz r oc r, while comoving temperature drops oc r . Since r cannot increase 
beyond r max ~ 77 ~ E /M c 2 , which is achieved at a radius rjr\ ~ 77, beyond this 
radius the flow begins to coast, with E ~ 77 ~ constant 90 ). However, the observed 
7-ray spectrum observed is generally a broken power law, i.e., highly nonthermal. 
The optically thick e ± 7 fireball cannot, by itself, produce such a spectrum (it would 
tend rather to produce a modified blackbody, 88 )' 39 )). In addition, the expansion 
would lead to a conversion of internal energy into kinetic energy of expansion, so 
even after the fireball becomes optically thin, it would be highly inefficient, most 
of the energy being in the kinetic energy of the associated protons, rather than in 
photons. 

The most likely way to achieve a nonthermal spectrum in an energetically effi- 
cient manner is if the kinetic energy of the flow is re-converted into random energy 
via shocks, after the flow has become optically thin m ) . Two different types of shocks 
may arise in this scenario. In the first case (a) the expanding fireball runs into an 
external medium (the ISM, or a pre-ejected stellar wind 111 )' 66 )' 48 )' 123 ). The second 
possibility (b) is that 112 )' 87 -*, even before external shocks occur, internal shocks de- 
velop in the relativistic wind itself, faster portions of the flow catching up with the 
slower portions. This is a completely generic model, which is independent of the 
specific nature of the progenitor, which has been successful in explaining the major 
observational properties of the gamma-ray emission, and is the main paradigm used 
for interpreting the GRB observations. 

External shocks will occur in an impulsive outflow of total energy E in an 
external medium of average particle density n Q at a radius 

r dec ~ (3E /4TTn m p c 2 V 2 ) 1/3 ~ lO 17 ^ 3 ^ 1 / 3 ^ 2 / 3 cm , (3-2) 

and on a timescale td ec ~ r dec/ {cE 2 ) ~ 3 x 10 2 E^ 3 n 1 ^f]2 ^ s > where rj = E = lO 2 ^ 
is the final bulk Lorentz factor of the ejecta. Variability on timescales shorter than 
tdec may occur on the cooling timescale or on the dynamic timescale for inhomo- 
geneities in the external medium, but generally this is not ideal for reproducing 
highly variable profiles 125 ). However, it can reproduce bursts with several peaks 94 ) 
and may therefore be applicable to the class of long, smooth bursts. 

The same behavior far with comoving temperature oc r _1 , followed by sat- 
uration r max ~ 77 at the same radius r/77 ~ 77 occurs in a wind scenario 89 ), if one 
assumes that a lab-frame luminosity L Q and mass outflow M Q are injected at r ~ ri 
and continuously maintained over a time t w \ here 77 = L /M c 2 . In such wind model, 
internal shocks will occur at a radius 112 ) 

Tdis ~ CtvarV 2 ~ 3 X l0 14 t var T]% CHI, (3-3) 

on a timescale t w 3> t var ~ r,ns/{crj 2 ) s, where shells of different energies Arj ~ 77 
initially separated by ct v (where t v < t w is the timescale of typical variations in 
the energy at n) catch up with each other. In order for internal shocks to occur 
above the wind photosphere r p h ~ M ax / {^m p cr 2 ) = 1.2 x 10 14 L5377^~ 3 cm, but 
also at radii greater than the saturation radius (so that most of the energy does not 
come out in the photospheric quasi-thermal radiation component) one needs to have 
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7.5 x lO 1 ^ 5 ^ 1 / 5 <??3 x 10 2 L^tva/ 4 ■ This type of models have the advantage 112 ) 
that they allow an arbitrarily complicated light curve, the shortest variation timescale 
t V ar ^10 -3 s being limited only by the dynamic timescale at 77, where the energy 
input may be expected to vary chaotically. Such internal shocks have been shown 
explicitly to reproduce (and be required by) some of the more complicated light 
curves 125 )' 54 )' 98 ) (see however 20 )). 

A potentially valuable diagnostic tool for the central engine of GRB is the power 
density spectrum (PDS). An analysis of a sample of 228 BATSE burst light curves 4 ) 
comes up with two interesting results, namely the logarithmic slope of the PDS be- 
tween 10~ 2 and 2 Hz is approximately -5/3, and there is a cutoff of the average PDS 
above 2 Hz. While these experimental results need further verification by indepen- 
dent data analysis groups, the experience with radio and accreting pulsars, black hole 
binaries and QPO sources indicates that time series analyses can provide valuable 
constraints on physical source mechanisms, e.g. 5 \ Using a simple kinematical model 

Fig. 2. Average power density spec- 
trum Pf of simulated bursts from 
internal shocks, compared with the 
observed PDS (thick line) deter- 
mined by 4 \ The average is done 
over 112 bursts with a square-sine 
modulated wind of random period 
between t w /4 and t w , with a power- 
law luminosity distribution and a 
constant comoving rate density evo- 
lution satisfyng the observed logN- 
logP constraints (Spada, Panaitescu 
& Meszaros 1999). 



for the shell ejection and collision and assuming the BATSE radiation to be due to 
synchrotron and/or IC scattered radiation, 100 )' 132 ) have calculated the light curves 
and PDS expected for a range of total burst energies and for a total mass ejected and 
bulk Lorentz factor distribution compatible with the internal shock scenario (Figure 
2). The distance distribution also contributes to the observed PDS, and this is taken 
into account by using a distribution compatible with the observed logN-logP counts. 
For optically thin winds, a slope approaching -5/3 requires a non-random Lorentz 
factor distribution, e.g. with an asymmetrical time modulation so as to produce 
a larger number of collisions at low frequencies. A cutoff at high frequencies (~ 2 
Hz) can be understood in terms of shocks which increasingly occur below the scat- 
tering photosphere of the outflow, or a deficit of energy in short pulses due to the 
modulation of the Lorentz factors favoring shocks arising further out. 

An interesting development is that significant fraction of bursts appear to have 
low energy spectral slopes steeper than 1/3 in energy 101 )' 18 ). This has motivated 
consideration of a thermal or nonthermal 59 -*' 60 ) comptonization mechanism. While 
an astrophysical model where this mechanism would arise naturally has been left 
largely unspecified, internal shocks which lead to pair formation 3S ) could have a 
self-regulated low pair temperature favoring comptonization. There is also growing 
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Fig. 3. Luminosity per decade xL x vs. x = hu/m e c 2 for two values of r) = L/Mc 2 and marginal 
(left) or large (right) pair compactness. T: thermal photosphere, PHC: photospheric comp- 
tonized component; S: shock synchrotron; C: shock pair dominated comptonized component 
(Meszaros & Rees, 1999b). 

evidence that the apparent clustering of the break energy of GRB spectra in the 
50-500 keV range may not be due to observational selection 101 )' 14 )' 21 \ Models to 
explain this, e.g., through a Compton attenuation model 15 ) require reprocessing by 
an external medium whose column density adjusts itself to a few g cm -2 . More 
recently a preferred break has been attributed to a blackbody peak at the comoving 
pair recombination temperature in the fireball photosphere 26 -*. In this case a steep 
low energy spectral slope is provided by the Rayleigh-Jeans part of the photosphere, 
and the high energy power law spectra and GeV emission require a separate explana- 
tion. In order for such photospheres to occur at the pair recombination temperature 
in the accelerating regime requires an extremely low baryon load. For very large 
baryon loads, a related explanation has been invoked 137 ), considerng scattering of 
photospheric photons off MHD turbulence in the coasting portion of the outflow, 
which upscatters the adiabatically cooled photons up to the observed break energy. 
These ideas have been synthesized 77 -* to produce a generic scenario in which the 
presence of a photospheric component as well as shocks subject to pair breakdown 
can produce steep low energy spectra and preferred breaks. 

§4. The Simple Standard Afterglow Model 

The dynamics of GRB and their afterglows can be understood in a fairly simple 
manner, independently of any uncertainties about the progenitor systems, using 
a generalization of the method used to model supernova remnants. The simplest 
hypothesis is that the afterglow is due to a relativistic expanding blast wave, which 
decelerates as time goes on 70 ) . The complex time structure of some bursts suggests 
that the central trigger may continue for up to 100 seconds, the 7-rays possibly 
being due to internal shocks. However, at much later times all memory of the initial 
time structure would be lost: essentially all that matters is how much energy and 
momentum has been injected; the injection can be regarded as instantaneous in 
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the context of the much longer afterglow. As pointed in the original fireball shock 
paper 111 ), the external shock bolometric luminosity builds up as L oc t 2 and decays 
as L oc t~( 1+q \ At the deceleration radius (3-2) the fireball energy and the bulk 
Lorentz factor decrease by a factor ~ 2 over a timescale t<i ec ~ r^ ec /(ci~' 2 ), and 
thereafter the bulk Lorentz factor decreases as a power law in radius. This is 

r oc r-n oc t-9/( 1 + 2 f) , r oc t x '^\ (4-1) 

with g is 3 for the radiative, 3/2 for the adiabatic regime (in which pr^T ~ constant 
or pr s r 2 ~ constant). 

The most obvious radiation mechanism is synchrotron, whose peak frequency 
in the observer frame is v m oc rB'^ 2 , and both the comoving field B' and electron 
Lorentz factor 7 are expected to be proportional to F 66 \ As r decreases, so will 
v m , and the radiation will move to longer wavelengths. Looking at the implications 
of this for the forward blast wave, 86 )' 49 ) discussed the possibility of detecting at late 
times a radio or optical afterglow of the GRB. A more detailed treatment of the 
fireball dynamics indicates that approximately equal amounts of energy are radiated 
by the forward blast wave, moving with ~ r into the surrounding medium, and by 
a reverse shock propagating with r r — 1 ~ 1 back into the ejecta 66 ). The electrons 
are shocked to much higher energies in the forward shock than in the reverse shock, 
producing a two-step synchrotron spectrum which during the deceleration time td ec 
peaks in the optical (reverse) and in the j/X (forward) 68 )> 69 ). Detailed calculations 
and predictions of the time evolution of such a forward and reverse shock afterglow 
model ( 70) ) preceded the observations of the first afterglow GRB970228 ( 17 )> 141 )), 
which was detected in 7-rays, X-rays and several optical bands, and was followed up 
for a number of months. 
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Fig. 4. The observed light 
curves of the afterglow 
of GRB 970228 at var- 
ious wavelenghts, com- 
pared 152 ' to the simple 
blast wave model predic- 
tions of 70) . 



The simplest spherical afterglow model concentrates on the properties of the 
forward blast wave only. The flux at a given frequency and the synchrotron peak 
frequency decay at a rate 70 - ) ' 76 -* 

F v OC t [3-2*(l-2/3)]/(l+2 fl ) ; ( 4 . 2 ) 
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where g is the exponent of r (equ. [4-1]) and (3 is the photon spectral energy 
slope. The decay rate of the forward shock F v in equ. (4-2) is typically slower than 
that of the reverse shock 70 \ and the reason why the "simplest" model was stripped 
down to its forward shock component only is that, for the first two years 1997- 
1998, afterglows were followed in more detail only after the several hours needed 
by Beppo-SAX to acquire accurate positions, by which time both reverse external 
shock and internal shock components are expected to have become unobservable. 
This simple standard model has been remarkably successful at explaining the gross 
features and light curves of GRB 970228, GRB 970508 (after 2 days; for early rise, 
see §5) e.g. 152 )' 135 )> 146 )' 114 ) (see Figure 4). 

This standard afterglow model produces at any given time a three-segment power 
law spectrum with two breaks. At low frequencies there is a steeply rising syn- 
chrotron self-absorbed spectrum up to a self-absorption break v a , followed by a 
+1/3 energy index spectrum up to the synchrotron break u rn corresponding to the 
minimum energy 7 m of the power-law accelerated electrons, and then a — (p — l)/2 
energy spectrum above this break, for electrons in the adiabatic regime (where 7 _p 
is the electron energy distribution above 7 m ). A fourth segment and a third break 
is expected at energies where the electron cooling time becomes short compared to 
the expansion time, with a spectral slope — p/2 above that. With this third "cool- 
ing" break uj,, first calculated in 74 ) and more explicitly detailed in 126 ), one has what 
has come to be called the simple "standard" model of GRB afterglows. One of the 
predictions of this model 70 ) is that the relation between the temporal decay index 
a, for g = 3/2 in r oc r~ g , is related to the photon spectral energy index (3 through 

F u oc Fv 13 , with a = (3/2)/? . (4-3) 

This relationship appears to be valid in many (although not all) cases, especially after 
the first few days, and is compatible with an electron spectral index p ~ 2.2 — 2.5 
which is typical of shock acceleration, e.g. 146 ^' 126 )' 151 \ etc. As the remnant expands 
the photon spectrum moves to lower frequencies, and the flux in a given band decays 
as a power law in time, whose index can change as breaks move through it. For the 
simple standard model, snapshot overall spectra have been deduced by extrapolating 
spectra at different wavebands and times using assumed simple time dependences 
147), 151) These can be used to derive rough fits for the different physical parameters 
of the burst and environment, e.g. the total energy E, the magnetic and electron- 
proton coupling parameters €b and e e and the external density n Q (see Figure 5). 

The simplest afterglow model is based on the following assumptions: a) A single 
value of E and r o = rj is used, b) the external medium n ext is homogeneous, c) 
the accelerated electron spectral index p, the magnetic field and electron to proton 
equipartion ratios eb and e e do not change in time, d) the expansion is relativistic 
and the dynamics are given by T oc r - 3 / 2 (adiabatic), d) the outflow is spherical, 
e) the observed radiation is characterized by the scaling relations along the line 
of sight. These assumptions, even if correct over some range, clearly would break 
down after some time. Estimates for the time needed to reach the non-relativistic 
expansion regime are typically < month(s) ( 143 - ) ), or less if there is an initial radiative 
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Fig. 5. Snapshot spectrum of 
GRB 970508 at t = 12 
days and standard after- 
glow model fit 151) . 



10 12 14 16 IB 
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regime r oc r -3 . However, even when electron radiative times are shorter than 
the expansion time, it is unclear whether a regime r oc r~ 3 should occur, since it 
would require strong electron-proton coupling 74 ) . As far as sphericity, the standard 
model can be straightforwardly generalized to the case where the energy is assumed 
to be channeled initially into a solid angle f2j < 4ir e7 \ In this case 117 )' 118 ) a 

change occurs after r drops below ^J 1 ^ 2 , after which the side of the jet becomes 
observable, and soon thereafter one expects a faster decay of r if the jet starts 
to expands sideways, leading to a decrease in the brightness. A calculation based 
on the sideways expansion, using the usual scaling laws for a single central line of 
sight 118 ) leads then to a steepening of the light curve. Until recently, no evidence 
for a steepening could be found in afterglows over several months. E.g., in GRB 
971214 110 ), a snapshot standard model fit and the lack of a break in the late light 
curve could be, in principle, interpreted as evidence for lack of a jet, leading to 
an (isotropic) energy estimate of 10 53 5 ergs. While such large energy outputs are 
possible in either NS-NS, NS-BH mergers 71 ) or in hypernova/collapsar models 91 )' 108 ) 
using MHD extraction of the spin energy of a disrupted torus and/or a central fast 
spinning BH, it is worth stressing that what these snapshot fits constrain is only 
the energy per solid angle 75 \ Also, the expectation of a break after some weeks or 

— 1/2 

months (e.g., due to r dropping either below a few, or below f2j ) is based upon 
the simple impulsive (angle-independent delta or top-hat function) energy input 
approximation. The latter is useful, but departures from it would be natural, and 
certainly not surprising. In fact, as discussed below, tentative evidence for beaming 
in one obejct has recently been reported 58 )' 33 )' 16 ), but it is difficulty to ascertain, 
and could be masked by a number of commonly expected effects. 

§5. More Realistic ("Post-standard") Afterglow Models 

Despite the success of the simple standard model, realistically one could expect 
any of several fairly natural departures from this to occur. One of these is that the 
emitting region, as seen by the observer, should resemble a ring 147 )' 95 )' 124 ). This 
effect may, in fact, be important in giving rise to the radio scintillation pattern seen in 
several afterglows, since this requires the emitting source to be of small dimensions, 
which is aided if the emission is ring-like, e.g. in the example of GRB 970508 148 ) . 
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Another likely complication is that the afterglows should show a diversity in their 
decay rates, not only due to different (3 but also from the possibility of a non- 
standard relation between the temporal decay index a and the spectral energy index 
f3, different from equ. (4-3). 

The most obvious departure from the simplest standard model occurs if the 
external medium is inhomogeneous: for instance, for n ex t °c r~ d , the energy con- 
servation condition is r 2 r 3 ~ d ~ constant, which changes significantly the temporal 
decay rates 74 ). Such a power law dependence is expected if the external medium 
is a wind, say from an evolved progenitor star as implied in the hypernova scenario 
(such winds are generally used to fit supernova remnant models). Another obvious 
non-standard effect, which it is reasonable to expect, is departures from a simple 
impulsive injection approximation (i.e. from a delta or top hat function with a sin- 
gle value for E Q and r o ). An example is if the mass and energy injected during the 
burst duration t w (say tens of seconds) obeys M(> r) oc r~ s , E(> r) oc i.e. 
more energy emitted with lower Lorentz factors at later times (but still shorter than 
the gamma-ray pulse duration). This would drastically change the temporal decay 
rate and extend the afterglow lifetime in the relativistic regime, providing a late 
"energy refreshment" to the blast wave on time scales comparable to the afterglow 
time scale 113 ). These two cases lead to a decay rate 

r -(3-d)/2 -n ext(xr -d. 

r -3/(2+s) .£ (> r) K r l- S W 

Expressions for the temporal decay index a(f3,s,d) in F u oc t a are given by 74 )' 113 ), 
which now depend also on s and/or al (and not just on (3 as in the simple standard 
relation of equ. (4-3). The result is that the decay can be flatter (or steeper, depending 
on s and d) than the simple standard a = (3/2)/3. A third non-standard effect, which 
is entirely natural, occurs when the energy and/or the bulk Lorentz factor injected 
are some function of the angle. A simple case is E Q oc 9~i , r o oc 9~ k within a range 
of angles; this leads to the outflow at different angles shocking at different radii 
and its radiation arriving at the observed at different delayed times, and it has a 
marked effect on the time dependence of the afterglow 74 \ with a = a((3,j, k) flatter 
or steeper than the standard value, depending on j, k. Thus in general, a temporal 
decay index which is a function of more than one parameter 

F v oat a v 13 ,with a = a(fi,d,s,j,k,- ■ ■) , (5-2) 

is not surprising; what is more remarkable is that, in many cases, the simple standard 
relation (4-3) is sufficient to describe the gross overall behavior at late times. 

Strong evidence for departures from the simple standard model is provided by, 
e.g., sharp rises or humps in the light curves followed by a renewed decay, as in GRB 
970508 ( 102 )> 106 )). Detailed time-dependent model fits") to the X-ray, optical and 
radio light curves of GRB 970228 and GRB 970508 show that, in order to explain 
the humps, a non-uniform injection or an anisotropic outflow is required. These 
fits indicate that the shock physics may be a function of the shock strength (e.g. 
the electron index p, injection fraction £ and/or e^, e e change in time), and also 
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indicate that dust absorption is needed to simultaneously fit the X-ray and optical 
fluxes. The effects of beaming (outflow within a limited range of solid angles) can 
be significant 97 ^, but are coupled with other effects, and a careful analysis is needed 
to disentangle them. 

One consequence of "post-standard" decay laws (e.g. from density inhomo- 
geneities, non-uniform injection or anisotropic outflow) is that the transition to a 
steeper jet regime r < 9J 1 ~ few can occur as late as six months to a year after 
the outburst, depending on details of the energy input. This transition is made 
more difficult to detect by the fact that, as numerical integration over angles of the 
ring-like emission 96 ) show, the transition is very gradual and the effects of sideways 
expansion effects are not so drastic as inferred 118 ) from the scaling laws along the 
central line of sight. This is because even though the flux from the head-on part 
of the remnant decreases faster, this is more than compensated by the increased 
emission measure from sweeping up external matter over a larger angle, and by the 
fact that the extra radiation, arising at larger angles, arrives later and re-fills the 
steeper light curve. The inference (e.g. 110 )' 118 )) that GRB 970508 and a few other 
bursts were isotropic due to the lack of an observable break is predicated entirely on 
the validity of the simplest standard fireball assumption. Since these assumptions 
are drastic simplifications, and physically plausible generalizations lead to different 
conclusions, one can interpret the results of 110 )' 118 ) as arguments indicating that 
post-standard features are, in fact, necessary in some objects. 

§6. Prompt multi-wavelength flashes, reverse shocks and jets 

Prompt optical, X-ray and GeV flashes from reverse and forward shocks, as well 
as from internal shocks, have been calculated in theoretical fireball shock models for 
a number of years 68 )' 69 )' 93 )' 70 )' 127 ), as have been jets (e.g. 64 )' 67 )' 69 ), and in more 
detail 117 )' 99 )' 97 )' 118 )). Thus, while in recent years they were not explicitly part of 
the "simple standard" model, they are not strictly "post" -standard either, since they 
generally use the "standard" assumptions, and they have a long history. However, 
observational evidence for these effects were largely lacking, until the detection of a 
prompt (within 22 s) optical flash from GRB 990123 with ROTSE by 2) , together with 
X-ray, optical and radio follow-ups citekul99,gal99,fru99,and99,cas99,hjo99. GRB 
990123 is so far unique not only for its prompt optical detection, but also by the fact 
that if it were emitting isotropically, based on its redshift z = 1.6 58 )' 1 ) its energy 
would be the largest of any GRB so far, 4 x 10 54 ergs. It is, however, also the 
first (tentative) case in which there is evidence for jet-like emission 58 )' 33 )' 16 ). An 
additional, uncommon feature is that a radio afterglow appeared after only one day, 
only to disappear the next 36 )' 58 ). 

The prompt optical light curve of GRB 990123 decays initially as oc t~ 2 - 5 to oc 
£-1.6 2)^ muc j 1 steeper than the typical oc t~ 1A of previous optical afterglows detected 
after several hours. However, after about 10 minutes its decay rate moderates, and 
appears to join smoothly onto a slower decay rate oc i~ L1 measured with large 
telescopes 36 )' 58 )' 33 )' 16 ) after hours and days. The prompt optical flash peaked at 9- 
th magnitude after 55 s 2 ) , and in fact a 9-th magnitude prompt flash with a steeper 
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decay rate had been predicted more than two years ago 70 ), from the synchrotron 
radiation of the reverse shock in GRB afterglows at cosmological redshifts (see also 
127 )). An optical flash contemporaneous with the 7-ray burst, coming from the 
reverse shock and with fluence corresponding to that magnitude, had also been 
predicted earlier 68 )' 69 ) . An origin of the optical prompt flash in internal shocks 70 )' 76 ) 
cannot be ruled out yet, but is less likely since the optical light curve and the 7-rays 
appear not to correlate well 128 )' 36 ) (but the early optical light curve has only three 
points). The subsequent slower decay agrees with the predictions for the forward 
component of the external shock 70 )' 128 )' 76 ). 

The evidence for a jet is possibly the most exciting, although must still be con- 
sidered tentative. It is based on an apparent steepening of the light curve after about 
three days 58 )' 33 )' 16 ). This is harder to establish than the decay of the two previous 
earlier portions of the light curve, since by this time the flux has decreased to a level 
where the detector noise and the light of the host galaxy become important. How- 
ever, after correcting for this, the r-band data appears to steepen significantly. (In 
the K-band, where the noise level is higher, a steepening is not obvious, but the issue 
should be settled with further Space Telescope observations). If real, this steepening 
is probably due to the transition between early relativistic expansion, when the light- 
cone is narrower than the jet opening, and the late expansion, when the light-cone 
has become wider than the jet, leading to a drop in the effective flux 117 )' 58 )' 76 )' 118 ). 
A rough estimate leads to a jet opening angle of 3-5 degrees, which would reduce the 
total energy requirements to about 4 x 10 52 ergs. This is about two order of magni- 
tude less than the binding energy of a few solar rest masses, which, even allowing for 
substantial inefficiencies, is compatible with currently favored scenarios (e.g. 108 )> 62 )) 
based on a stellar collapse or a compact binary merger. 

§7. Location and Environmental Effects 

The location of the afterglow relative to the host galaxy center can provide clues 
both for the nature of the progenitor and for the external density encountered by the 
fireball. A hypernova model would be expected to occur inside a galaxy in a high 
density environment n a > 10 3 — 10 5 cm~ 3 . Most of the detected and well identified 
afterglows are inside the projected image of the host galaxy 9 ), and some also show 
evidence for a dense medium at least in front of the afterglow ( 85 ) ) . For a number of 
bursts there are constraints from the lack of a detectable, even faint, host galaxy 13 °) , 
but at least for Beppo-SAX bursts (which is sensitive only to long bursts tb>20 s) 
the success rate in finding candidate hosts is high. 

In NS-NS mergers one would expect a BH plus debris torus system and roughly 
the same total energy as in a hypernova model, but the mean distance traveled 
from birth is of order several Kpc 12 ), leading to a burst presumably in a less dense 
environment. The fits of 151 ) to the observational data on GRB 970508 and GRB 
971214 in fact suggest external densities in the range of n Q = 0.04-0.4 cm~ 3 , which 
would be more typical of a tenuous interstellar medium. These could be within the 
volume of the galaxy, but for NS-NS on average one would expect as many GRB 
inside as outside. This is based on an estimate of the mean NS-NS merger time of 10 8 
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years; other estimated merger times (e.g. 10 7 years, 140 )) would give a burst much 
closer to the birth site. BH-NS mergers would also occur in timescales ,<10 7 years, 
and would be expected to give bursts well inside the host galaxy ( 12 ) ; see however 34 ) ) . 
In at least one "snapshot" standard afterglow spectral fit for GRB 980329 116 ) the 
deduced external density is n Q ~ 10 3 cm -3 . In some of the other detected afterglows 
there is other evidence for a relatively dense gaseous environments, as suggested, 
e.g. by evidence for dust 115 ) in GRB970508, the absence of an optical afterglow and 
presence of strong soft X-ray absorption 40 )' 79 ) in GRB 970828, the lack an an optical 
afterglow in the (radio-detected) afterglow ( 136 )) of GRB980329, and spectral fits 
to the low energy portion of the X-ray afterglow of several bursts 85 ). The latter 
observations may be suggestive of hypernova models 91 )' 34 ), involving the collapse of 
a massive star or its merger with a compact companion. 

One important caveat is that all afterglows found so far are based on Beppo-SAX 
positions, which is sensitive only to long bursts if>,>20 s 46 ). This is significant, since 
it appears likely that NS-NS mergers lead 62 ) to short bursts with tb<W s. To make 
sure that a population of short GRB afterglows is not being missed will probably 
need to await results from HETE 42 ) and from the planned Swift 134 ) mission, which 
is designed to accurately locate 300 GRB/yr. 

An interesting case is the apparent coincidence of GRB 980425 with the unusual 
SN Ib/Ic 1998bw 35 ), which may represent a new class of SN 47 )' 10 ) . If true, this could 
imply that some or perhaps all GRB could be associated with SN Ib/Ic 145 ), differring 
only in their viewing angles relative to a very narrow jet. Alternatively, the GRB 
could be (e.g. > ) a new subclass of GRB with lower energy £ 7 ~ 10 48 (i7 i /4vr) erg, 
only rarely observable, while the great majority of the observed GRB would have 
the energies ~ 10 54 ( fij/Air) ergs as inferred from high redshift observations. The 
difficulties are that it would require extreme collimations by factors 10~ 3 — 10~ 4 , and 
the statistical association is so far not significant 51 ) . However, two more GRB light 
curves may have been affected by an anomalous SNR (see, e.g. the review of 150 )). 

The environment in which a GRB occurs should also influence the nature of the 
afterglows in other ways. The blast wave and reverse shock that give rise to the 

— 1/3 

X-rays, optical, etc occur over timescales proportional to td ec <x n ex t ( ec L u -[3 - 2]) 
which is longer in lower density environments, so for the same energy the flux is 

1/2 

lower, roughly F v oc E n ext , contributing also to make afterglows in the intergalactic 
medium harder to detect. However, in addition to affecting broad-band fluxes, one 
may also expect specific spectral signatures from the external medium imprinted in 
the X-ray and optical continuum, such as atomic edges and lines 6 )' 104 )' 73 ). These 
may be used both to diagnose the chemical abundances and the ionization state (or 
local separation from the burst), as well as serving as potential alternative redshift 
indicators. (In addition, the outflowing ejecta itself may also contribute blue-shifted 
edge and line features, especially if metal-rich blobs or filaments are entrained in 
the flow from the disrupted progenitor debris 72 ), which could serve as diagnostic 
for the progenitor composition and outflow Lorentz factor). To distinguish between 
progenitors (§2), an interesting prediction ( 73 ); see also 37 )' 13 )) is that the presence of 
a measurable Fe K-a X-ray emission line could be a diagnostic of a hypernova, since 
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Fig. 6. Spectrum (top) of a 
hypernova funnel model 
for various observer times, 
showing (bottom) the to- 
tal and Fe light curves 
and equivalent widths 
(Weth, Meszaros , Kali- 
man & Rees 1999), for 
with R=1.5x 10 16 cm, 
n = 10 10 cm~ 3 , and Fe 
abundance 10 2 times so- 
lar. 
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in this case one may expect a massive envelope at a radius comparable to a light- 
day where vr<l, capable of reprocessing the X-ray continuum by recombination 
and fluorescence. Detailed radiative transfer calculations have been performed to 
simulate the time-dependent X/UV line spectra of massive progenitor (hypernova) 
remnants 149 ), see Figure 6. Two groups 107 )' 156 ) have in fact recently reported the 
possible detection of Fe emission lines in GRB 970508 and GRB 970828. 



§8. Conclusions 



Gamma Ray bursts provide perhaps the most direct method for studying strong 
gravity phenomena near black holes, and with over 4,000 bursts detected, they would 
represent the most abundant supply of high quality photon data on black holes. The 
fireball shock model of gamma-ray bursts has proved quite robust in providing a 
consistent overall interpretation of the major features of these objects at various 
frequencies and over timescales ranging from the short initial burst to afterglows ex- 
tending over many months. The standard internal shock scenario is able to reproduce 
the properties of the 7-ray light curves, while external shocks involving a forward 
blast wave and a reverse shock are successful in reproducing the afterglows observed 
in X-rays, optical and radio. The "simple standard model" of afterglows, involving 
four spectral slopes and three breaks is quite useful in understanding the 'snapshot' 
multiwavelength spectra of most afterglows. However, the effects associated with 
a jet-like outflow and the possible differential beaming at various energies requires 
further investigations, both theoretical and observational. Caution is required in 
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interpreting the observations on the basis of the simple standard model. For in- 
stance, more detailed numerical models, as opposed to the more common analytical 
scaling law models, show that the contributions of radiation from different angles 
and the gradual transition between different dynamical and radiative regimes lead to 
a considerable rounding-off of the spectral shoulders and light-curve slope changes. 
Time-dependent multiwavelength fits ") of some bursts also indicate that the param- 
eters characterizing the shock physics change with time, and a non-standard relation 
between the spectral and temporal decay slope observed in several objects, e.g. GRB 
990123 5S ), may provide evidence 76 ) for "post-standard" effects in such bursts. 

Significant progress has been made in understanding how gamma-rays can arise 
in fireballs produced by brief events depositing a large amount of energy in a small 
volume, and in deriving the generic properties of the long wavelength afterglows that 
follow from this. There still remain a number of mysteries, especially concerning the 
identity of their progenitors, the nature of the triggering mechanism, the transport 
of the energy, the time scales involved, and the nature and effects of beaming. How- 
ever, even if we do not yet understand the details of the gamma-ray burst central 
engine, it is clear that these phenomena are among the most powerful transients in 
the Universe, and they could serve as powerful beacons for probing the high redshift 
(z > 5) universe. The modeling of the burst mechanism itself, as well as the resulting 
outflows and radiation, will continue to be a formidable challenge to theorists and to 
computational techniques. Nonetheless, the collective theoretical and observational 
understanding is vigorously advancing, and with dedicated new and planned obser- 
vational missions under way, further significant progress may be expected in the near 
future. 
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